Density-functional theory of inhomogeneous electron systems in thin quantum wires 
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Motivated by current interest in strongly correlated quasi-one-dimensional (ID) Luttinger liquids 
subject to axial confinement, we present a novel density-functional study of few-electron systems 
confined by power-low external potentials inside a short portion of a thin quantum wire. The 
theory employs the ID homogeneous Coulomb liquid as the reference system for a Kohn-Sham 
treatment and transfers the Luttinger ground-state correlations to the inhomogeneous electron sys- 
tem by means of a suitable local-density approximation (LDA) to the exchange-correlation energy 
functional. We show that such ID-adapted LDA is appropriate for fluid-like states at weak cou- 
pling, but fails to account for the transition to a "Wigner molecules" regime of electron localization 
as observed in thin quantum wires at very strong coupling. A detailed analyzes is given for the 
two-electron problem under axial harmonic confinement. 

PACS numbers: 71.10.Pm,71.10.Hf,71.15.Mb 
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I. INTRODUCTION 

One-dimensional (ID) quantum many-body systems 
of interacting fermions have attracted theoretical and 
experimental interest for more than fifty years^. Con- 
trary to what happens in higher dimensionality, these 
systems cannot be described by the conventional Landau 
theory of normal Fermi liquids 2 due to a subtle inter- 
play between topology and interactions. The appropriate 
paradigm for ID interacting fermions is instead provided 
by the Luttinger-liquid concept introduced by Haldane 
in the early eighties 3 . 

Strongly correlated ID systems that are nowadays 
available for experiment range from ultra-cold atomic 
gases^ to electrons in single-wall carbon nanotubes^ and 
in semiconductor quantum wires&i. Chiral Luttinger 
liquids at fractional quantum-Hall edges 8 - also provide 
an example of ID electronic conductors and have been 
the subject of intense experimental and theoretical stud- 
i cc& i P . In many experimental situations the translational 
invariance of the fluid is broken by the presence of in- 
homogeneous external fields. Examples are the confin- 
ing potential provided by magnetic and optical traps for 
ultra-cold gases^ and the barriers at the end of a quan- 
tum wire segment in cleaved edge overgrowth samples 7 . 
These strong perturbations induce the appearance of a 
new length scale and can cause novel physical behav- 
iors relative to the corresponding unperturbed, Galileian- 
invariant model system. 

Of special relevance to the present work are the stud- 
ies carried out in Ref. 0, where momentum-resolved tun- 
neling experiments between two closely situated paral- 
lel quantum wires have been carried out to probe the 
phenomenon of spin-charge separation in a Luttinger liq- 
uid— In these experiments a top gate is used to deplete 
the central portion of one of the two wires, thus locally 
decreasing the electron density, and a dramatic transition 
is observed when the electron density is reduced below a 
critical value. There is strong evidence that in this regime 
the electrons in the depleted wire segment are separated 



by barriers from the rest of the wirei 2 -, and it is sug- 
gested that the electrons in the segment are localized by 
the combined effect of the barriers and of the electron- 
electron interactions. The magnetic-field dependence of 
the tunneling conductance for a field perpendicular to the 
plane of the wires provides a direct probe of the many- 
body wavefunction of the localized electron s 12 ' 13 , offering 
the possibility to investigate systematically the role of in- 
teractions in creating exotic phases of matter in reduced 
dimensionality. In fact, the experimental parameters in 
Ref. lD are such that the electrons in the wire segment 
are in the strong-coupling regime. An exact diagonaliza- 
tion study was carried out in Ref. [12| for a number N of 
electrons up to 4. 

A powerful theoretical tool to study the interplay 
between interactions and inhomogeneity from external 
fields of arbitrary shape is density-functional theory 
(DFT), based on the Hohenberg-Kohn theorem and the 
Kohn-Sham mapping 14 . Many-body effects enter DFT 
via the exchange-correlation (xc) functional, which is of- 
ten treated by a local-density approximation (LDA) re- 
quiring as input the xc energy of a homogeneous refer- 
ence fluid. For 2D and 3D electronic systems the un- 
derlying reference fluid usually is the homogeneous elec- 
tron liquid (EL) , whose xc energy is known to a high de- 
gree of numerical precision from quantum Monte Carlo 
(QMC) studies 1 ^. Several density-functional schemes 
have also been proposed for strongly correlated ID sys- 
tern s 16 ' 17 ' 18 ' 19 ' 20 and in the case of the ID Luttinger 
liquid with repulsive contact interactions, where the xc 
energy of the homogeneous fluid is exactly known from 
Bethe-Ansatz solutions, tests of the LDA have been car- 
ried out against QMC data—. 

In the present work we test a novel LDA xc functional 
to treat few-electron systems confined by power-law po- 
tentials inside a segment of a quantum wire. The ho- 
mogeneous reference system that we adopt is the ID 
EL with Coulomb interactions, previously studied by 
a number of author o 21 ' 22 ' 23 and very recently evalu- 
ated by a novel lattice-regularized Diffusion Monte Carlo 
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method^. The correlation energy determined in this lat- 
ter study is used in our LDA calculations, whereas earlier 
DFT-based studies of ID "quantum dots'— have used 
the correlation energy of a 2D EL. While such choice 
can be justified for thick wires, a II? reference fluid is 
more appropriate to treat inhomogeneous electron sys- 
tems in ultrathin wires of our present interest (see also 
the discussion given in Ref. 23). We nevertheless find 
that the ID-adopted LDA is unable to describe the tran- 
sition of the confined electrons from a fluid-like state to 
the localized "Wigner-like" state that is observed to oc- 
cur as the coupling strength is increased. In essence, 
the confining potential pins the phase of density oscil- 
lations in much the same way as an impurity inserted 
into the infinitely extended II? fluid docs in producing 
Friedel oscillations in the surrounding electron density. 
However, a cross-over from a 2k-p to a 4fcp periodicity 
occurs in these oscillations with increasing coupling in 
the Luttinger liquid. We proceed in the later part of the 
paper to give a detailed analysis of this transition in the 
case of two electrons subject to axial harmonic confine- 
ment in a wire segment. From previous work on the two- 
particle problem with contact repulsive interactions^ we 
presume that a local spin-density approximation could 
help in transcending the limitations of the LDA. 

The outline of the paper is briefly as follows. In Sect. UT1 
we introduce the Hamiltonian that we use for the system 
of present interest, and in Sect. IIIII we describe our self- 
consistent DFT approach and the LDA that we employ 
for the xc potential. In Sect. lIVl we report and discuss our 
main results for the fluid state at weak coupling, while in 
Sect. [V] we focus on the two-electron problem. Finally, 
Sect. IvTI summarizes our main conclusions. 



II. THE MODEL 

We consider N electrons of band mass m confined in- 
side an axially symmetric quantum wire. The transverse 
confinement is provided by a tight harmonic potential 
with angular frequency uix, 



Vx(x,y) = ^muj 2 ± (x 2 +y 2 ) . 



(1) 



The electrons are also subject to a longitudinal poten- 
tial V ext (z) along the wire axis. In the ID limit (see 
below) the transverse motion can be taken as frozen 
into the ground state of the 2D oscillator, ip(r±) = 
(27rfe 2 )- 1 / 2 exp[-ri/(46 2 )] with b 2 = h/(2mujx). The 
parameter b thus measures the transverse wire radius. 
On integrating out the transverse degrees of freedom one 
ends up with the effective ID Hamiltonian 



n 



where 



2m 



2^ 9 ? + 9 2^ 



Vb(\zi-Zj\) + ^2 V e xt(Zi) , (2) 



is the renormalized interelectron potential^ 2 -. Here k is 
a background dielectric constant and erfc(x) is the com- 
plementary error functional. It is easy to check that the 
potential in Eq. ((3]) becomes purely Coulombic at large 
distance 28 , i>b(z) — > e 2 /(n\z\) for \z\ — > oo. At zero inter- 
electron separation the electron-electron potential goes to 
a positive constant. Equation (|3|) yields a linear approach 
to a constant, the cusp being an artifact of wavefunction 
factorization 2 ^. 

The last term in Eq. gives the coupling of the 
electrons to the axial external potential and, following 
Tserkovnyak et ai.— , we consider power-law potentials of 
the type 



V cxt (z) = Vp\zf 



(4) 

). For (3 = 2 
raw 2 z 2 1 2 with 



with (5 > 2 and V = 2 l3+1 h 2 /(mL 2+l3 
the confinement is harmonic, V cxt (z) = 

angular frequency wn = Ah/ (mL 2 ), while V ex t(z) becomes 
a square well of size L in the limit /3 — > +oo. 

Choosing L/2 as the unit of length and 2H 2 / (mL 2 ) as 
the unit of energy, the Hamiltonian becomes 



(5) 



with_x = 2z/L, A = _ ^/irL 2 /(Ab)_ and T(x) = 
exp[L 2 £ 2 /(16& 2 )]erfc[La;/(46)]. Here L = L/a B and 
b = b/a-Qi clb — h 2 K,/(me 2 ) being the effective Bohr ra- 
dius. We see from Eq. (O that the physical properties 
of the system are determined by the four dimensionless 
parameters N, b, f3, and L. Note that while T(x) is con- 
trolled only by the ratio L/b, the parameter A contains 
two powers of L and one power of b. Electron-electron in- 
teractions are expected to become dominant in ultrathin 
wires with b < 1 and for weak confinements (L 1). In 
the experiments of Ref. (with as — 9.8 nm for GaAs) 
b w 1 and L w 100, so that the electrons are in a strong- 
coupling regime (A « 4 x 10 3 ). In fact, the electron- 
electron coupling is also influenced by the exponent /3, 
which determines the spill-out of the electron density and 
hence the system diluteness. For given b and L, harder 
boundaries (larger j3) imply a more efficient confinement 
(i.e. higher average density) and thus reduce the role of 
the many-body interactions. 



III. DENSITY-FUNCTIONAL APPROACH 

Within the Kohn-Sham version of DFT the ground- 
state density riGs(z) is calculated by self-consistently 
solving the Kohn-Sham equations for single-particle or- 
bitals ip a (z), 



v b (z) 



2 Kb 



exp[z 2 /(46 2 )]erfc[z/(2&)] 



(3) 



-L^ +VKs[nGs]{z) 



Ifa(z) = £ a ip a (z) (6) 
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with Vks(z) = vn(z) + v xc (z) + V C xt(z) , together with the 
closure 



n G s(z) 



r Q \<p a {z)Y 



(7) 



Here the sum runs over the occupied orbitals and the 
degeneracy factors T a satisfy the sum rule J2 a ^ a = N. 
The hrst term in the effective Kohn-Sham potential is 
the Hartree term 



+ OC 



dz'vb{\z - z'\)n GS (z') , (8) 



while the second term is the xc potential, defined as the 
functional derivative of the xc energy E xc [n] evaluated at 
the ground-state density profile, i> xc = SE xc [n]/ Sn(z)\ G s- 
The total ground-state energy of the system is given by 



dzv xc [n GS }(z)n G s(z) 
- - / dz dz'v b (\z - z'\)n GS (z)n GS (z') 



E xc [n G s] ■ 



(9) 



Equations j6|) and ((7]) provide a formally exact scheme 
to calculate n G s (z) and E G g , but E xc and i> xc need to be 
approximated. 

As mentioned above in Sect. HI in this work we have 
chosen the ID EL, described by the Hamiltonian ([2]) with 
V cx t{z) — 0, as the homogeneous reference fluid. In the 
thermodynamic limit and in the absence of spin polar- 
ization this model is described by two dimcnsionless pa- 
rameters only, r s and b. Here r s — (2naB) 1 is the usual 



Wigner-Seitz dimensionless parameter, defined in terms 
of the average ID density n. We adopt the LDA func- 
tional 



+ OO 



dzn(z)e*r(r s (z)) (10) 



with r s {z) = [2n GS {z)a B }- 1 and £ xc om (^) = 4° m (r s ) + 
eJ? om (r s ). The exchange energy e x om of the ID EL (per 
particle) is calculated from 



-horn / 



+ 00 



p-Mi)i s o(q) 

Zir 



(11) 



where Vb(q) — (e 2 / k) exp (q 2 b 2 )F<i(q 2 b 2 ) is the Fourier 
transform of the interaction potential, with Ei(a;) 
being the exponential integral^!, and So(q) is the 
noninteracting-gas structure factor (So(q) — q/(2ky) for 
q < 2kp and 1 elsewhere). The correlation energy eJ? om 
determined by Casula et al£± is given by the parametriza- 
tion formula 



_hom 



(rs) = 



A + Br] + Cr 2 



In (1 + Dr s +Er]). 



(12) 

in units of the effective Rydberg e 2 /(2naB). The val- 
ues of the seven parameters in this expression are re- 
ported in Table IV of Ref. [24] for several values of b 
in the range 0.1 < b < 4. As discussed in Ref. 0, 
Eq. (fi"2")) incorporates the exactly-known weak-coupling 
limit (r s — > 0) and fits very well their numerical data in 
the range 0.05 < r s < 50. Finally, the LDA xc potential 
is calculated from Eq. (fTU]) as 



.LDA 



n G s}(z) 



5n 



GS 



(13) 



r„^[2?iGs(z)aB] _1 



We have calculated numerically the derivative of the ex- 
change energy as 



1 f 1 

= -7T2 — / d( i {q~ l ) v b(q) 

+ rV/ 1 «»- 1 ^- (14) 

2r s a B Jq or s 



where q — q/(2kp). Notice that Vb(q) is rvdependent. 



IV. NUMERICAL RESULTS FOR THE FLUID 
STATE 

We have solved numerically the self-consistent scheme 
given by Eqs. ©-([I]) using the LDA xc potential in 
Eq. (|13p . Our main numerical results for the density 
profile n G s(z) of even numbers of electrons in a weak- 
coupling regime are summarized in Figs. [TJ3J 

In the homogeneous limit the ID hypothesis (i.e. a sin- 
gle transverse subband occupied) requires that the Fermi 
energy £p = h 2 kp/(2m), with fcp = nn/2 — 7r/(4r s aB), 
be smaller than the transverse energy Hluj_. This trans- 
lates into the inequality r s > 7T&/4, involving r s and the 
wire radius b in units of the Bohr radius. In our calcu- 
lations we have checked that the minimum r s (z) defined 
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FIG. 1: (Color online) Top panel: Density profile ugs(z) (in 
units of 2/L) as a function of 2z/L for = 4, 6, and 8 elec- 
trons confined by an external potential with (3 = 2 and L = ab 
in a thin wire of radius b = 0.1o,b. Bottom panel: Same as in 
the top panel but for (3 = 6. 



FIG. 2: Top panel: Density profile ugs(z) (in units of 2/L) 
as a function of 2z/L for N = 6 electrons confined by an 
external potential with (3 = 2 and L/as = 1, 2 and 3 in a 
thin wire of radius b = O.Iob. Bottom panel: Same as in the 
top panel but for (3 = 6 and L/cib = 2, 4 and 6. Results for 
the noninteracting system are also shown in both panels for 
comparison. 



by the local density tigs(z) satisfies the ID hypothesis 
for each set (iV, b, (3, L) of parameters. 

In Fig. Q]we report the density profiles for N — 4, 6, 
and 8 electrons in the case of a thin wire with radius 
b = 0.1 as and a confinement with L — ae, correspond- 
ing to A ~ 4. We see from this figure that for these 
system parameters the ground state is fluid-like with 
N/2 distinct maxima, corresponding to Friedel-like oscil- 
lations with wave number 2kp where the effective Fermi 
wavenumber = irn/2 is determined by the average 
density n in the bulk of the trap. In Fig. [5] we show 
the evolution of the density profile with increasing L 
for N — 6 electrons confined in a thin wire of radius 
b = 0.1 ob, and in Fig. [3] we show the evolution of the 
density profile with increasing b for fixed L — 2a^- The 
role of electron-electron interactions becomes more im- 
portant with increasing L or decreasing b (for L = 6ae 
and b = OTae for example, we have A ~ 160) and leads to 
a decrease in the amplitude of the Friedel-like oscillations 
and to a broadening of the density profile. 

In Fig. 2] we report the dependence of the ground- 
state energy Eqs and of the stiffness d 2 Eos/dN 2 — 
[E GS {N + 2) + E GS (N - 2) - 2E GS (N)}/4: on the elec- 
tron number TV, for different types of confining potential. 
The behavior of these quantities is easily understood in 
the noninteracting case. In harmonic confinement the 
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FIG. 3: (Color online) Density profile riGs(z) (in units of 
2/L) as a function of 2z/L for iV = 6 electrons confined by an 
external potential with (3 = 2 and L = 2as, for two values of 
the wire radius. Results for the noninteracting system have 
also been included for comparison. 
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FIG. 4: (Color online) Top panel: Ground-state energy [per 
particle and in units of 2fi 2 /(mL 2 )] as a function of the elec- 
tron number N > 2 for various values of /3 and 6 = O.larj. 
Filled symbols correspond to the interacting case, empty sym- 
bols to noninteracting case. For — 2 we have chosen L — as, 
while for /3 = +oo L — 5o,b. Bottom panel: Thermodynamic 
stiffness d 2 EGs/dN 2 [in units of 2fi 2 /(mL 2 )] as a function of 
N for the same system parameters as in the top panel. The 
lines are just guides for the eye. 



single-particle spectrum is given by = tkJu(i + 1/2) 
with i = 0,1,2,... and thus the ground-state energy is 



i = huj\\N 2 /4, implying a constant 
= ftwii/2. In the case (3 = +oo, 



stiffness d 2 E GS /dN 2 

instead, e< = h 2 n 2 i 2 j '(2mL 2 ) with i = 1,2,3, ... and thus 
Egs(N) = h 2 TT 2 N(N + 1)(N + 2)/(24mL 2 ), implying 
a linear stiffness d 2 E GS /dN 2 = h 2 n 2 (N + l)/(4mi 2 ). 
We are instead unable to calculate the addition energy 2 ^ 
(chemical potential) /i = Eqs(N) — Eqs(N — 1), as it 
requires knowledge of the ground-state energy for sys- 
tems having odd numbers of electrons and hence a finite 
spin polarization. The spin-polarization dependence of 
the correlation energy of the ID EL is presently not yet 
available. 

Whereas the above results refer to a fluid-like weak- 
coupling regime, one should expect real-space quasi- 
ordering to set in at strong coupling, and this should 
be signaled by the so-called "2fcp — * 4kp crossover" in 
the wave number of Friedel oscillations. This cross-over 
is not predicted by the LDA xc functional in Eq. (fl"3)) . 
In Section [V] we study in detail this crossover for N = 2 
harmonically-trapped electrons, a problem which is eas- 
ily solvable numerically to any desired degree of accuracy 




2z Ic i/ L 

FIG. 5: (Color online) The effective potential V(z re \) [in units 
of 2fi 2 /(ml 2 )] as a function of 2z Te \/L for b = O.Iob- 



(see also the work of Szafran et a/.—) 



V. THE TWO-PARTICLE PROBLEM AND THE 
FAILURE OF THE LDA AT STRONG COUPLING 

After a canonical transformation to centre-of-mass and 
relative coordinates and momenta [Z = [z\ + Z2)/2,P = 
Pi +P2 and Zi e i = zi — Z2,p = (pi — P2}/2], the Hamilto- 
nian for two harmonically trapped electrons in a thin wire 
can be written as 7i = Hcm(Z, P) + 7^ re i(z re i,p). Here, 
the centre-of-mass Hamiltonian Hcm — P 2 /(2M) + 
Moj 2 Z 2 /2 describes a ID harmonic oscillator of mass 
M = 2m, while the relative-motion Hamiltonian H. re \ = 
p 2 /m + V(Zrel) describes a particle of mass m/2 in the 
potential V(z ro i) = mu 2 z 2 cl l '4 + Vb(z rc \). This potential 
is plotted in Fig. [5] for two values of the trap frequency 
W|| = 4/i/(m£ 2 ). 

In the spin-singlet case the spatial part of the ground- 
state wavcfunction is written as 



*gs(zi, z 2 ) = A/"exp (-Z 2 /a 2 ,)(p ro \(z Ic i) , 



(15) 



where Af is a normalization constant, ay = y/h/(mu>u), 
and ip re i(z re i) is the symmetric ground-state wavefunction 
for the relative- motion problem with energy e r , which 
can be numerically found by solving the single-particle 
Schrodinger equation 



h 2 d 2 

m dz 2 cl 



+ V(z lel ) 



¥>rclOrcl) = £r¥>rel(Zrel) • (16) 



An illustration of I^qsC 2 : z ')\ 2 for two values of L is re- 
ported in Fig. [6] The "molecular" nature of the ground 
state is evident at strong coupling. 

The ground-state density profile can be found from 



dz'\* GS (z,z')\ 2 , 



(17) 
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where the normalization constant J\f is chosen according 
to J^^dzriQs(z) = 2. Numerical results are shown in 
Fig. [JJ (left panels) in comparison with the LDA profiles. 
Note that the double-peak structure in I^gsC^, z')\ 2 a t 
weak coupling (left panel in Fig. [6]) is lost in the corre- 
sponding ground-state density. While at weak coupling 
(L = as) the agreement between the exact result and the 
LDA prediction is very satisfactory, at strong coupling 
(L = 4ob) the LDA is unable to reproduce the formation 
of a deep Coulomb hole yielding a density profile with a 
broad maximum at the trap center. 

One can also directly compare the LDA xc potential in 
Eq. (fLl)) with the exact one, which can be calculated from 
the exact density profile 3 -^ as summarized below. In the 
two-particle case there is only one Kohn-Sham orbital 
fKs(z) = \J riQs(z)/2, which satisfies the Kohn-Sham 
equation 



In Fig. [7] (right panels) we show a comparison be- 
tween the LDA xc potential in Eq. (flU)) , as obtained at 
the end of the Kohn-Sham self-consistent procedure, and 
the exact xc potential calculated from Eq. (f2"TJ)) with the 
ground-state density from Eq. (|17[) . Several remarks are 
in order here. As it commonly happens, the LDA po- 
tential has the wrong long-distance behavior: it decays 
exponentially because the density does so, while the ex- 
act xc potential decays like l/\z\. Nevertheless, at weak 
coupling the difference between the two potentials is well 
approximated by the constant 5, in the region where the 
density profile is different from zero, and this explains the 
satisfactory agreement between the exact and the LDA 
profiles. It is finally evident how in the strong-coupling 
regime the LDA potential is instead very different from 
the exact one. 



2m dz 2 



+ V K s[riGs](z) 



tPKs(z) = £ks<Pks(z) . (18) 



Solving this equation for v xc we find 

t u \ h2 d 2 Lp KS (z) 
v X c[nGs\(z) = £ks + ~ TT T~? Vcxt(z) 



2mifiKs(z) dz 2 



- v H (z). 
or, more explicitly, 

v xc [nGs]{z) = £ks + 



d 2 ^/n G s(z) 



+ 00 



2my / n GS (z) dz 2 
dz' v b (\z - z'\)n GS (z') . 



(19) 



V cxt {z) 



(20) 



The exact Kohn-Sham eigenvalue £ks can be proven to 
be equal to the energy e r of the relative motion. The 
approximate Kohn-Sham eigenvalue, instead, differs from 
e r : for example, for L — ob we find S = e^g A — e r ~ 
0.46 [2h 2 /(mL 2 )}. 



An xc functional that embodies the 2kp — > 4fcp 
crossover and is capable of describing inhomogeneous 
Luttinger systems at strong repulsive coupling is thus 
required. In Ref. [26| we have proposed a simple xc func- 
tional which is able to capture the tendency to antiferro- 
magnetic spin ordering. The idea consists in two steps: 
(i) one adds an infinitesimal spin-symmetry-breaking 
field to the Hamiltonian; and (ii) one resorts to a lo- 
cal spin-density approximation (LSDA) within the frame- 
work of spin-density functional theory. Earlier exact di- 
agonalization and configuration- interaction studies of ID 
quantum dot o 30 ' 32 have shown that, while for even num- 
ber of electrons the local spin polarization is everywhere 
zero in the dot, one can still observe antifcrromagnetic 
correlations at strong coupling by looking at the spin- 
resolved pair correlation functions. This suggests that an 
LSDA approach may indeed prove useful at strong cou- 
pling. Unfortunately, a knowledge of the ground-state 
energy of the homogeneous II? EL in the situations with 
iV T ^ N t is still lacking. 
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FIG. 7: Left panels: Density profile t%gs(z) (in units of 2/L) as a function of 2z/L for N — 2 electrons confined by a harmonic 
potential with L = as (top) and L = 4as (bottom) in a thin wire of radius b = 0.1o,b. The exact results (filled circles) are 
compared with the LDA results (solid line). Right panels: The LDA xc potential from Eq. (|13|l (solid line) is compared with 
the exact xc potential calculated from Eq. (|20p (filled circles). 



VI. CONCLUSIONS 

In summary, we have carried out a novel density- 
functional study of a few isolated electrons at zero net 
spin, confined by power-law external potentials inside a 
short portion of a thin semiconductor quantum wire. The 
theory employs the quasi-one-dimensional homogeneous 
electron liquid as the reference system and transfers its 
ground-state correlations to the confined inhomogeneous 
system through a local-density approximation to the ex- 
change and correlation energy functional. 

The local-density approximation gives good-quality re- 
sults for the density profile in the liquid-like states of 
the system at weak coupling, a precise test against ex- 
act results having been presented in the case of N = 2 
electrons. However, it fails to describe the emergence of 
electron localization into Wigner molecules at strong cou- 



pling. The fact that strong-coupling antifcrromagnetic 
correlations are hidden in the inner-coordinates degrees 
of freedom, as suggested by Szafran et al£° , indicates 
that a local spin-density approximation, or even non-local 
functionals based on the spin-resolved pair correlation 
functions^, are needed. The class of density-functional 
schemes for "strictly correlated" electronic systems re- 
cently proposed by Perdew et alM* may also be useful in 
treating the Wigner-molecule regime. 
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